Entanglement and Spontaneous Symmetry Breaking in Quantum Spin Models 
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It is shown that spontaneous symmetry breaking does not modify the ground-state entanglement 
of two spins, as defined by the concurrence, in the XXZ- and the transverse field Ising-chain. Cor- 
relation function inequalities, valid in any dimensions for these models, are presented outlining the 
regimes where entanglement is unaffected by spontaneous symmetry breaking. 
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Entanglement is a property of a quantum state shared 
between two or more parties. It is denned with the aim 
of capturing the essential quantum non-locality encoded 
in the state. While for a long time interests in entan- 
glement stemmed from the opportunity to understand 
fundamental concepts in quantum mechanics such as the 
EPR-paradox and violation of Bell-inequalities, recent in- 
terest in entanglement comes from its use as a resource 
for performing tasks not possible by classical means. 

From an applied viewpoint it is thus worth quantify- 
ing the degree of entanglement in natural systems, such 
as solid-state materials. A number of studies have been 
devoted to quantifying such "natural" entanglement in 
states of simple models describing idealized quantum 
magnets, such as the XXZ- and transverse field Ising- 
models pi y, S U III L| • However the possible alteration of 
entanglement by spontaneous symmetry breaking (SSB) 
were not discussed in any of these works, although the 
need for such a study was mentioned in Ref. [f| . SSB hap- 
pens invariably in real materials described by a Hamilto- 
nian possessing a global symmetry, thus it is important to 
investigate whether or not the entanglement calculated in 
the symmetric ground state(s) is changed by SSB. Here 
we find that entanglement is not modified by SSB for the 
XXZ- and transverse field Ising-chains. 

It is a fundamental requirement of entanglement that 
it cannot on average be created by mixing classically two 
quantum states. Classical mixing is incoherent, and so 
one should not gain more entanglement from this mix- 
ing than what is already encoded in the states that are 
being mixed. Thus if we consider the density matrix 
p = (p+ + p~)/2 put together by an equal mixture of the 
broken symmetry states, p± (not necessarily pure), the 
entanglement E(p) must satisfy 



E{p) < \ (E(p, 



E(p-)) = E(p±) 



(1) 



where we have assumed in the last equality that the en- 
tanglement in the two different symmetry-broken ground 
states are equal, the states being related by a global 
change of basis states. Thus the entanglement in the 
broken state cannot be smaller than in the symmetric 
state. 



Here we focus on the entanglement of two specific spins 
(or qubits) in the ground-state of a quantum system. The 
state of two spins i and j in the ground state of a quantum 
system is described in terms of the reduced density ma- 
trix pij obtained by tracing over all spins in the ground 
state except the two spins i and j. Writing it out ex- 
plicitly in the standard basis {| TT), I Tl), I IT), I II)} onc 
has 



({P}P}) (P}o-) iarpj) (a-aT)\ 



Pij 



(phi) 

?i pi 



(2) 



where P T = ±(1 + a z ), P l = ±(1 - a z ) and a ± = 
t- (a x ± ia v ). The brackets denote ground state expec- 
tation values and a are the Pauli matrices. 

Recently entanglement between blocks of consecu- 
tive spins have been studied using the von Neumann 
entropy0 which is the conventional measure of bipar- 
tite entanglement in pure states. There it is essential that 
the ground-state is pure as the von Neumann entropy is a 
concave function violating Eq. JQ| . This requirement be- 
comes very restrictive when considering the state of just 
two spins tracing out the rest. Then even though the 
ground-state of the full system is pure, pij is generally 
not. Thus a definition of entanglement in mixed states 
is needed. A reasonable definition is the entanglement of 
formation, Ef which is obtained by decomposing the 
mixed ensemble into pure states, and summing the en- 
tanglement in each of these pure states weighted by their 
probabilities of occurrence in the mixed ensemble. As 
there are many possible decompositions into pure states, 
Ef is taken to be the minimal value gotten by trying all 
possible decompositions. This definition of entanglement 
clearly satisfies Eq. Q as the minimal decompositions 
of p + and p_ into pure states also is a decomposition of 
p, but not necessarily the minimal one. Furthermore Ef 
in pure states reduces to the conventional definition of 
entanglement. 

The minimization procedure over all possible decom- 
positions is however quite difficult to handle practically. 
It is therefore of great value that there exist a closed- form 
expression for this minimum value for bipartite two-level 
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systems. This expression is in terms of the concurrence 
C , which is defined in terms of the spectrum of the ma- 
trix PijPij where p = a\ ® a v . p* a\ ® a v - , the time-reversed 
density matrix. The complex conjugation refers to the 
standard basis used in Eq. J2J. Let A/ be the eigenvalues 
of pp so that Ai > A 2 > A3 > A4. Then the concurrence 
C is 



C = max 



{», 



Ai — v A2 — v A3 — 



A 4 } (3) 



and the entanglement of formation is Ef = — xlog 2 x — 
(1 — x) log 2 (l — x), where x = 1/2 + vl — C 2 /2. Because 
Ef is a monotonous function of C with Ef(C = 0) = 
and Ef(C = 1) = 1 we will hereafter for simplicity discuss 
the concurrence instead of the entanglement of formation. 
In Ref. 10] it was shown that the concurrence as defined 
above is the largest convex function reducing to the pure 
state concurrence, thus C satisfies Eq. Q . 

Before discussing mixed ensembles of general states we 
will first discuss the ensemble obtained by an equal clas- 
sical mixture of two pure states with equal concurrences. 
That is p = {p + + /Q-)/2 and p± = |a±)(a±|. The two 
states \a±) can be thought of as different realizations 
of a broken global symmetry. Denoting the concurrence 
c = |(a_|_|o!_|_)| = |(a_|S_)|, and the overlap (a+\a-) = d 
the square roots of the eigenvalues of pp are |c± |<2||/2. 
Thus 



C(p) = mm{c, \d\} . 



(4) 



So the concurrence of the symmetric mixed state equals 
the concurrence of the broken symmetry state only when 
|(a+|d+)| < |(a + |aT_)|. From this it becomes clear that 
the question of the influence of SSB on concurrence can- 
not be answered in general terms, one needs to look at the 
specifics of the states involved. We will therefore inves- 
tigate two specific models; the XXZ- and the transverse 
field Ising-model. 

The Hamiltonian for the XXZ-model is 



E{- 

(id) 



Atrfo^} , (5) 



where the sum is taken over all nearest neighbor sites on 
a lattice which for simplicity is taken to be bipartite. The 
Hamiltonian is real and invariant under a [/(l)-rotation 
about the spin z-axis. This symmetry being continuous 
can only be broken in dimensions higher than one and for 
|A| < 1. A global 7T-rotation about the spin x (or y) axis 
is also a symmetry (Z2). This symmetry can be broken 
in one dimensional systems when |A| > 1. Let us first 
consider the breaking of the ^-symmetry assuming the 
f7(l)-symmetry is unbroken. The reduced density matrix 
then reads 



Pij 



/AO \ 

B C 

C G 

\ D J 



(6) 




FIG. 1: The concurrence as function of A for different separa- 
tions between sites i and j for a 128 site XXZ-chain with(filled 
symbols, only for A > 1) and without (open symbols) a sym- 
metry breaking staggered field h z = 0.001. Fhe inverse tem- 
perature is 1/T = 80. For A > 1 the open and filled symbols 
lie on top of each other. 



where the symbols A . . . G expressed in terms of correla- 
tion functions can be read off Eq. ©. With a matrix of 
this form the concurrence is C — 2 • max{0, \C\ — \J AD} 
which in terms of correlation functions is 



C 



- • max { 0, (7) 

\{afa*) + (afaj)! - ^ {I + (a.f a*))* - (af + } 



The breaking of the Zi symmetry manifests itself as 
(cr z ) 7^ 0. The values of the other correlators are un- 
changed by SSB as erferj, rf<jj, vfo-j are all invariant 
under the Z2 transformation. 

For A < — 1 SSB does not affect the concurrence as it is 
zero in both cases, the ground state being either a prod- 
uct state or a sum of product states. Thus the interesting 
region is the antiferromagnetic Ising region A > 1 where 
(erf) = ±m, the sign depending on which sublattice spin 
i is located on. For i and j on different sublattices , the 
concurrence is not affected by the symmetry-breaking as 
then (erf + <j z a ) = 0. For i and j on the same sublattice 
knowledge about the values of the correlation functions 
is needed. To obtain these we have performed Monte 
Carlo simulations using the Stochastic Series Expansion 
technique [ll| with directed-loop updates^]. The cor- 
relation functions were obtained directly from the loop 
movement in space-time in analogous manner to that de- 
scribed in Ref. |13| . In order to simulate the broken sym- 
metry phase a small symmetry breaking staggered field 
were imposed. The results can be seen in Fig. ^ which 
shows the concurrence as a function of A for different 
separations between i and j both with and without a 
symmetry-breaking staggered field along the spin z-axis. 
As can be seen the symmetry-breaking field has no effect 
on the concurrence. For i and j nearest neighbors the 
concurrence reaches a maximum of 0.386 at A = 1 in 
agreement with Ref. [lj], while for bigger distances be- 
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tween i and j the concurrence drops rapidly to zero long 
before A reaches 1. 

In dimensions greater than one the [/(l)-symmetry is 
broken when |A| < 1. If the ^-symmetry is obeyed, the 
density matrix takes the form 



pij 



( A a a f \ 

a B C a 

a C B a 

\f a a A J 



(8) 



where a = (of)/4 = (oJ)/4 and / = (of of - ofo?)/4 
are the only values affected by the breaking of the U(l)- 
symmetry. (The case when a magnetic field is present 
along the spin-z axis is discussed at the very end of this 
article.) The square roots of the eigenvalues of pp written 
in terms of the correlation functions are in this case 



u ± = 7 



V± 



^(l + ^foj)) 2 - 4(^)2 ± |(ofoJ) 
±((oM) + (ofof))|. 



\ i j ) 



(9) 



In the symmetric state (a a 



and (of cr? 



Ma]) 



The condition that p^ has non-negative eigenvalues im- 
plies that whenever the concurrence is non-zero in the 
symmetric state, u + is the largest of the square roots of 
the eigenvalues. Now assume that the symmetry is bro- 
ken in a continuous manner and that the concurrence is 
non-zero in the symmetric state, then it + remains the 
largest until another eigenvalue takes over that role, at 
which point the concurrence becomes zero. However, if 
this were to hap pen the fundamental criterion Eq. 
would be violated[la|- Thus we conclude that u+ remains 
the largest even when the symmetry is broken. 

The concurrence in the symmetry broken state then 
takes the t/(l)-invariant form 



C = \ ({of of, 



-i 



(10) 



iff <^> + (of(7?) > (ofo?) - 1 and (o?**) > (ofof). 
If these inequalities are not satisfied, the concurrence 
will take on different ^7(l)-7io?i-invariant expressions in 
the broken state. However, of these only IjlUI) reduces 
smoothly (it is constant) to the non-zero value gotten 
at a, f — > 0, thus we conclude that the expression ijTTJI) is 
always valid in a regime close to the symmetric state pro- 
vided the concurrence in the symmetric state is non-zero. 
The results for the case of antiferromagnetic interactions 
in the spin-xy-plane can be obtained by changing signs on 
<j x and <j v on one sublattice in the ferromagnetic results 
described here. 

The Hamiltonian for the transverse field Ising model is 



H 



E 



o, a a + 



(11) 
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FIG. 2: Ising chain with 256 sites in a transverse field with and 
without a symmetry breaking field. The inverse temperature 
is 1/T = 80. The top panel show the nearest-neighbor concur- 
rence C calculated from the measured nearest neighbor den- 
sity matrix vs. A for different values of the symmetry breaking 
field. The open symbols means h x = while the closed sym- 
bols is for h x = 0.001. The inset in the top panel shows a 
blow-up of the region close to the critical point A = lambda c 
as well as data for h x = 0.01 (filled triangles). The bottom 
panel show the transverse and longitudinal magnetization for 
the case h x = 0.001 as functions of A. The dotted lines are ex- 
act results at T — for zero symmetry breaking field|l6|. the 
value of (a x ) is gotten from (crfcrj) x (cr x ) 2 for \i — j\ — > oo. 



This model has less symmetries than the XXZ-model. In 
fact the only symmetry is a global 7r-rotation about the 
spin z-axis. This symmetry being discrete is broken for 
A = l/(2h z ) > \ c in any dimensions. To gain insight we 
have performed Monte Carlo simulations both with and 
without a symmetry breaking field h x along the spin-x 
axis, see Fig. [21 The results without a symmetry break- 
ing field is in good agreement with the results obtained 
in Refs. 0,0- From the upper panel in Fig. [3 it is seen 
that the effect of the symmetry breaking magnetic field 
is to slightly lower the concurrence in the critical region 
around A c . This would seem to contradict Eq. l[T]l. how- 
ever from the inset in the upper panel we see that this 
effect diminishes as h x — ► 0, thus we interpret the low- 
ering as an effect of the small but finite h x , which will 
change the symmetry broken states themselves whenever 
the gap to excitations is small as is the case close to 
the critical point. The deviation from the ideal situation 
when h x just acts as a symmetry breaking field can also 
be seen in the lower panel where the magnetization is 
seen to deviate slightly from the exact result close to A c . 
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We conclude that the Monte Carlo results indicate that 
SSB does not alter the nearest-neighbor concurrence in 
the transverse field Ising-chain. 

To explain this in more details consider the density 
matrix for two spins in the transverse field Ising model 



(A 



Pij 



a a F\ 

a B C b 

a C B b 

\F b b 



(12) 



where the symmetry-breaking manifests itself in non- 
zero values of a and b reflecting the non-zero values 
of (<t ± ) and (ct ± ct z ). For the symmetric state where 
a = b = 0, the square root of the eigenvalues of pp 
are B ± C and V AD ± F, thus the concurrence is 

C = 2 • max jo, |C| - VAD, \F\ - s|. For general values 
of a, b the eigenvalues of pp is gotten by solving the de- 
tcrminantal condition which amounts to solving a cubic 
polynomial equation as the eigenvalue \B — C\ 2 factorizes 
out. Instead of solving for the eigenvalues we formulate a 
condition the coefficients go,gi and gi of this cubic poly- 
nomial must satisfy in order for the concurrence to stay 
constant. Let x 2 , y 2 and z 2 be the roots of the cubic poly- 
nomial in increasing order. Then the equations relating 
these eigenvalues to the coefficients are 



.90 = (xyz) , 

2 2 i 2 2 i 2 2 

9i = x y + x z + y z , 
92 



(13) 



2 i 2 i 2 

x + y + z . 



The condition that the concurrence stays constant is 



z — [x + y) = k, 



(14) 



where k is a constant equal to the a = b = value. 
Squaring Eq. I|14l) . rearranging, and squaring again one 
finds that the roots can be completely eliminated in favor 
of the coefficients, thus 



I 2 

2^^/50 = j (k 2 - 32) -gi- 



(15) 



If this equation holds for k constant as the coefficients g 
are varied, the concurrence is unaffected by SSB. Written 
out 

.go = [a 2 - &j6]/3 - Apva - Ap 2 5 + Av 2 ^ 

gi = a 2 + 2a(3 - Apv - 4 7 <5 (16) 

.92 = 2a + 13 

where a = F 2 + AD - 2ab, (3 = (B + C) 2 - 4a6, 7 = 
DF - 6 2 , 8 = AF - a 2 , p = aD - b(B + C — F) and 
v = a(B + C — F) — bA. Squaring Eq. fT5|l and inserting 
the expressions for the g's one can check that Eq. i|15|) 
holds for nonzero a and b iff V 'AD + F is the biggest of 



the square roots of the eigenvalues; i.e. k = 2F— (B + C). 
In terms of correlation functions this condition is 



-4(<7. 



z\2 



- 1 > 2 (o? <T 



y _2A 



(17) 



We have checked numerically that this inequality is sat- 
isfied for non-zero values of h z for the nearest neigh- 
bor correlators in the transverse field Ising chain with 
no symmetry-breaking field, thus SSB does not alter 
nearest-neighbor concurrence in this case. 

The reduced density matrix for the XXZ-model in the 
presence of an external magnetic field along the spin z- 
axis takes the same form as Eq. (1121 with the difference 
that F = f is not invariant under the ?7(l)-symmetry 
rotation about the spin z-axis. In the ?7(l)-symmetric 
case, a = b = f = and so whenever concurrence is 
non-zero |i?-(-C| 2 is the largest eigenvalue, implying that 
k = B+C— 2^/ AD. However with this value of k Eq. l(T5|) 
is not satisfied for non-zero values of a, b and /. Thus we 
conclude that in this case SSB is generally accompanied 
by a change in the concurrence. 

While we have given conditions for when the concur- 
rence is not changed by SSB in some specific models it 
would be nice if more general statements can be made. It 
is reasonable to expect that these will involve statements 
about the system's behavior under time-reversal. 

The author acknowledges helpful discussions with Alan 
Luther. 
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